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Method* and apparatus for measuring flissne histology 



This invention relates to methods end devices for fbe'non iixvasiver measurement of 
epithelial tissue histology and is parti oularly, but not exclusively, concerned with 
measuring sHn histology, These cpfttafol tissues include fee respiratory tract, Che genital 
tract, the gastrointestinal tract and the retina of the eye. 

The distribution of chromophobes within these tissues and flheir lesions is considered to 
be potentially ussM information for the screening of patients to identify those whs 
should be referred to an appropriate clinician for farther diagnosis through tissue biopsy 
o* other conventional tests- Devices to detect these distortions with non invasive 
techniques are therefore beneficial. Further, the simplification of such devices And the 
automation of the diagnostic process will make such devices available to wider 
populations. 

The present invention improves the measurement technique described in GB 9624003, i 
and related patents, and PCIYGB00/Q2124. The devices described in this application are 
improvements of those described in PCT/GBOl/01986, PCT/GBOl/0303 1 and 
GBO 112501.2. 

Figure 12 shows a cposs section of skin. The major structural layers are well faxown s 
comiSed layer 45, epidermis 46, dermo-epidernml jimction 47, papiftaay dermis 48, 
re^cute dermis 49, In normal skin melanocytes 50, the ceils responsible for the 
production of the pigment melanin, He in the epidermis. In some benign skin lesions, 
such as the compound nevus melanocytes can exist in a characteristic pattern in both the 
epidermis and the papillary dermis. In melanoma the melanocytes can be asymmetrically 
distributed in both the epidermis and papillary dennis. 

The present invention provides additional information on the distribution of the pigment 
melanin in the papillary dermis, and includes a method for analysing this data 
automatically. This method is also applicable to other epithelial tissues with similar 
layered structures with differing optical properties- 
Method to analyse presence and depth of chromophores in tissue 

The equation described in OB application 9624003,1 page 10 lines 9 & 10 defines the 
expected colouration of human skin exhibiting the dermal descent of melanin as; 



(Equation Z) 

Hie terms of this equation are defined in GB application 9624003.1 pages 6 to 10* The 
equation may be modified to suit other tissues and chromophores and in the following 



d u *tl, i .tt lt .0 m d m )m 
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description the terms 'skhi' and 'melanin* represent a particular example application of 
this method* The technique is particularly applicable to the ohomophore haemoglobin in 
the skin where this method can be readily applied. 

In the example of measurement of the location of melanin in the sldn this equation is 
plotted in figure 1 with two primaries PI & F2 chosen such that blood concentration is 
constant and melanin can vary boflh in concentration ajjd location in the skin. 

Point 1 is the position of infinite melanin concentration in any position in the skin, point 
2 is the point of no melanin in the skin. The line 5 connects points of increasing melanin 
concentration in the epidermis 45. 

Figure 12 shows that that there are three variables of imerest that define the distribution 
of melanin in the skccu The concentration of melanin in the epidermis, the concentration 
of melanin in the papillary dermis and the location of melanin in the papillary dermis. 
Two simplified cases are considered, the first where the concentration of melanin in the 
papillary dermis is assumed to be uniform, and secondly that the depth of dermal melanin 
is either zero or the depth of the psjpillary dermis itself. 

Standard dermal melanin concentration case 

In figure 1 point 3 i$ a point of chosen melanin concentration in the epidermis. As 
dermal melanin depth is increased (at standard concentration) the line 6 is generated. 
Point 4 represents the ma x imum possible depth that is the depth of the papillary dermis 
itself. This depth may have been standardised by the date calibratipn steps described in 
GB 9624003.1. Figure 3 shows the effect of repeating lime 6 with different amounts of 
dexmal melanin. Lines 5.10,1 U 12 are then plotted joining points on these lines where the 
same dermal melanin depth is used. Lines 5,10,1 1>12 define aress 7,8,9 that represent 
areas where colours of shallow, medium and deep dermal melanin depth will lie. ' 

Standard dermal melanin depth case 

j 

Figure 4 show? the same axes, curve 5 with end points 1 & 2 as described above. In this 
case a point of epidermal melanin concentration is chosen 13 and, two lines of low and 
high dermal melanin concentration are plotted with increasing depth 1 7, 16. The points 
at which the mTO.imum dermal melanin depth is reached (papillary dermis depth) on these 
lines are 14 and 15. 

Figure 5 shows these points 14, 15 joined to form lines of constant dermal melanin 
concentration at maximum depth 18, 19»20. These tines bound areas 21,22^3 that 
represent areas where colours of low, medium and high dermal melanin concentration 
will lie, assuming dermal melanin exists and it*s maximum possible depth. 

It can be semi from this analysis that a single point on these two axes is insufficient to 
folly define the dermal melanin depth, dermal melanin concentration and epidermal 
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S2S conc f n I tmt, ' on - The following invention overcomes this problem and allow* a 
clinically usefiil measurement to be made. 

S^S™ 1 J dU ^ ° f ^ Iailin COnCentratfon m fte epidermis, melanin concentration in 
the papllary demus and dermal melanin depth that typically occur in human skin and in 

™^5F W ^ 0bS T:f ^ ^ 21 ' 22 ' 23 » fi ** e ^ and me areas 7,8,9 in figure 3 
were similar even though they were plotted using different simplification as™tiolT 
These area, represent similar values of a m(dermal melanin depth, dermal ineS 
concentration, epidermal melanin conceniratito). This function can also be interpreted as 

^H^ m f OT J°A mEaSU ^ Qf ,faal meknin vo W - This apprtncJS^af be 
mtuabvly understood by considering figure 5 and figure 3 . Line 5 represents tissue with 
no dermal melanin. Areas 21 end 9 represent tissuewith either deep Jow^nSZfef 
melamn, or shallow high concentration melanin (approximating to tissue with a low 
dennal meiamn volume'). Whereas areas 23 and 7 represent tissue with deep, Hen 
vcwTr^ C«PP«>xlmating to tissue with high ^nnaTmelanin 

S^! } Jr?* O^ 8 ^^ spectral response of tissue where significant melanin 
has penetrated the papillary dermis. 'Dermal melanin volume' is clinically significant 
because m a typical melanoma lesion the melanin producing cells, melanocytes, are 
similar because they are all produced by reproduction from a single mutated celL and 
therefore changes in melanin concatenation are likely to be small compared to changes in 
the looahon of these cells. Therefore in such lesions areas with increased 'dermal melanin 
volume area areas with increased dermal melanin depth. Dermal melanin depth is a 
well-known diagnostic and prognostic indicator in melanoma. If 'dennal melanin 
ZOUS ^ ?Jf f^yfP*™* **oss an area of skin, the measurements can be 

2f n ^ a ^played as a visualisation approximating to a map of dermal 
melanin depth across the lesion. This can assists the clinician is deteanining the degree 
of chaos in the distribution of melanocytes within the lesion. 

A particular technique to make measurements of &(dernMil melanin depth, dermal 
melanin concentration, epidermal melanin concentration) or 'dermal melanin volume' is 
as toUows. A spectral measurement is made using a devioe as described in figures 8 9 & 
iw*™* 3 dsscrib^in PCT/GBOl/01986 orPCT/GBOl/03011, aridmemWements 
plotted on a graph indicatong the areas of similar values of fh(Dermal melanin depth, 
Dermal melanin concentration, Epidemiol melanin concentration) (figure 6) This 
measurement can men be categorised by fh(Dermal melanin depth, Dennal melanin 
concentration, Epidermal meianin concentration). 

A second particular technique to make measurements of fitfdermal melanin depth, dermal 
melanin concentration, epidermal melanin concentration) or 'dermal melanin volume' is 
as loiiows. A spectral measurement is made using a device as described in figures 8,9 & 
10 or devices described in PCT/GB01/0I986 orPCTVGBOl/0301 1 and plotted on the 
graph shown m figure 5 , This measurement is compared to lines of changing dermal 
meiamn concentration with constant dennal melanin depth. The line that intersects the 
measurement point provides a value of dennal melanin concentration which can be used 
as an approximation to the value of fitfdermal melanin depth, dermal melanin 
concentration, epidermal melanin concentration). Figure 13 describes this process. 
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la one particular implementation of the process described in figure 13 these 
measurements are taken over a series of adjacent 15 micron square pixels ™aH) mm 
area, and the spectral measurements are made are me percentage remittance of light ftom 
the tissue at three separate primary wavebands with peaks at 450nm 550wu and 650nm, 
each with a SOnm bandwim. 

Alternatively the measured point could be plotted on figure 3 where me measurement is 
compared to lines of changing dermal melanin depth wito constant dermal melanin 
concentration. The line that intersects the point provides a value of dermal melamn depth 
which can be used as an approximation to the vajue ofm(dermal melamn depth, dermal 
melanin concentration, epidermal melanin concentration). In this case porat 3 would 
represent a zero value, and point 4 a maximum value. 

If an array of spectral measurements are made and a corresponding array of fii(Dermal 
melanin depth, Dermal melanin concentration, Epidermal melanin concentration) is 
calculated, then a colour can be assigned to values of m(dermal melamn depth, dermal 
melanin concentration, epidermal melanin concentration) and the data represented as a 
false colour image. 

En a further invention me approximate value of fn(dermaj melanin depth, dermal melamn 
concentcation, epidermal melamn concentration) obtained ftom the above methods can be 
further analysed to enable suspicious lesion to be differentiated from benign ones. The 
following analysis technique can also be applied to arrays of other measurements made 
spatially across a tissue sample, 3^pariicmarmemeasDreroentsmaybe ofbloodor 
collagen concentration. 

In me particular case of skin lesiona, malignant sldn lesions often contain Hie pigment 
melanin in the papillary dermis of the skin. A device that detects the presence ofo^rmal 
melanin is effective at detecting maKgnant lesions and tnwh devices are described m our 
patent applications GB0016690.0 andGBOl 12501.2. However tteae devices suffer from 
false positives generated by benign lesions which also contain dermal melamn. It is weB 
known that several types ofbenign sldn lesions contain dermal melanin, mclumng toe 
compound nevus and blue nevus. Therefore to increase the specificity of a device to 
detect malignant lesions a method of differentiating benign lesions containing dermal 
melanin for malignant lesions is beneficial. 

A well known property of malignant lesions is toe asymmetric and ch ^l™**** 
lesion. This invention is a method of automatically measuring the chaotic growth of toe 
lesion. 

Figure 7 shows a diagram of a typical skin lesion viewed ftom above. The skin 26 
contains a roalanooytic lesion 25, containing epidermal melanin as well as a region 
containing both epidermal and dermal melanin 27. 
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This mveation creates two data arrays to describe the lesion derived from spectral 
measurements using the techniques described in the te*t above and our previous patents 
which are also detailed above. 

fix this particular embodiment one array contains values proportional to the total amount 
of melanin in the skin, and a second with values appmxfm&tibtig to fh(dennal melanin 
depth, dermal melanin concentration, epidermal melairin concentration). Whore no 
dermal melanin wrists this second array has zero values* 

A well -known property of skin lesions ia they have characteristic patterns in the 
distribution of melanin- Benign lesions have a variety of patieans* but importantly these 
patterns are similar in both tibte epiderrai^ and papillary dennis as melanocytes slowly 
descend from the dermal/epidermal junction into the papillary dermis. These melanocytes 
can also cease melanin production limiting dermal melanin to a shallow oven layer in the 
papillary dermis. Malignant lesions characteristically grow chaotically in three 
dimensions so the patterns of melanin distribution are significantly different in the 
epidermis and papillary dermis. This invention automatically detects these differences 
through a mathematical process- 
Figure 1 1 describes this process. In a particular implementatioxi it operates on two arrays 
of data generated from spectral measurements taken over an area of skin. One array Is of 
measurements of total melanin calculated using the method described in GB9624003 . 1 
and PCT/GB00/Q2124, and the second is an array of values of &(dedmal melanin depth, 
dermal melanin concentration, epidermal melanin concentration) or "dermal melanin 
volume* calculated from the process described in figure 13. 

In step 1 of figure 11 the lesion is identified by applying a. threshold to values of total 
melanin TM(x,y). to step 2 the textuie of the total melanin array within the lesion is 
measured. In one implementation the standard deviation SD is used, but other well 
known statistical measures could also be used such as range or fractal dimension. Bi step 
3 areas inside the lesion where dermal melanin is present are identified as those 
coordinates where the 'dermal melanin volume* DMV(x,y) is not zero. In step 4 the same 
texture measurement is made of DMVfoy). In this implementation the standard 
deviation SD is used again, but other well known statistical measures could also be used 
such as range or fractal dimension. In step 5 a ratio of the two texture measurements is 
made and in step 6 this is threshold is. applied to determine malignancy. 

An improvement to this sequence is to add a preliminary step which detects elements in 
the array where hairs ate present on the skin and marks them to be excluded fiom the 
calculation process. Erroneous data from the presence of hair has been found to disrupt 
this algorithm and provide unreliable results. 

A particular embodiment of a device to implement the above techniques and those 
described in GB9624003, FCT/GBOD/02124 and GBQ1 12501-2 is described below and in 
figures 8 and 9. This device consists of a light tight housing 30 containing a transparent 
optical window 29 which is placed in contact with the akin 28. Light of a specific 
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wavelength band is generated by light emitting diodes mounted on a ring 33. This ring 
may contain a number of groups of LEDs each with it's own emission waveband. In one 
particular embodiment there are four groups of LEDs, each group containing between 8 
and 24 individual LED dies, this number depei*ds on the brightness of the individual 
devices, where relatively dim devices are used in larger numbers to approximately 
equalise the illumination intensity across wave bands* In oneparticular embodiment the 
wavebands of each group of LEDs span between 400nm and 1 OOOnm, In one particular 
embodiment an additional reflector i$ mounted beside each LED to ensure a majority of 
the emitted light is reflected onto the window 19* hi one particular embodiment an 
additional filter is added over one or more of the LED groups to eliminate unwanted 
emissions in the mfta-red region of the spectrum, Ibis xfeg 33 is covered by a polarising 
filter 32. An important feature of the design of the LED system is that the LEDs are 
evenly spaced relative to the window 29 to ensure even iUujxrination of the skin 28. The 
LED illuminator is controlled by an electronic driver system 38 which receives control 
signals from the processing mejms 39. The electronic driver system responds to signals 
received from the processing means and switches the LED groups on and off in sequence 
as required to synchomise with the camera system 35 and 37.; In one embodiment the 
driver system can deliver different currents to each LED group, these currents can be 
generated according to data passed from the processing means 39 to Ihe Hit) control 
system 38. hi one embodiment the driver system can deliver diflterent currents to each 
LED to compensate for the differences in efficiency of each device. Light remitted from 
the skin 28 passes through a second polarisation filter 31 mounted such that its 
polarisation axis is at 90 degrees relative to the polarisation axis of Slter 32. This light 
then passes into the lens 34, and is focussed onto, a CCD $ensor array 35 . hi a particular 
embodiment the lens 34 is designed to focus affectively over, a wide waveband 
corresponding to the wavebands emitted by the USD illwdnator 33. In a particular 
embodiment this waveband is from 400 nm to 1 00Q m The CCD sensor is controlled 
by an electronic system 37 which converts the light intensity on each array element into a 
digital pixel value in an image array which is passed to the processing means 39. In a 
particular implementation the processing means can control the exposure time and gain of 
the CCD array 35 > The complete device is controlled by a control programme 
implemented by the processing means 39, This programme provides for a sequence of 
images to be acquired, each with a different illumination spectrum provided by the LED 
light source 33 and a corresponding exposure time- In on* particular embodiment eight 
individual images arc recorded, one with each of four wavelength bands with exposures 
set suph that normal skin can be recorded* On the three wavelength bands which lie in the 
visible spectrum an additional oyer exposed image is recorded with 4 times the standard 
exposure time for that wavelength band. The eighth image is taken with no illumination 
to determin the black values for each picture element. Tbfc image acquisition cycle is 
triggered by switch 36- The image data is then processed using algorithms described 
above. The process data is displayed on a screen 40. 

One particular embodiment of this device is illustrated in figure 9, in which items 29 to 
38 are contained within a handset 41 and the processing means and display means are^ 
provided by a laptop computer 42. This provides a compact and portable implemenation 
of this device, 
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In a second embodiment, illustrated in figure 10, all items 29 - 40 axe integrated into a 
hand held battery operated device. In this embodiment the output is simplified using the 
method described above and in figure 1 1 and an output is presented * $flia!1 LCD screen 
43. The optical window 29 is shown as 44* 
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Figure 3 
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Figure 10 
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Step 1 -Identify lesion boumJry 

Lesioxafoy) = ttue if TM(x,y) > threshold 
value 



Step 2 ~ Measure texture of total melanin 
distribution (TIM) in the lesion 

TTM « SD ofTMfoy) where 
lesionC^y^true 

. i 

Step 3 - Identify areas of dsraoal melanin i» 
lesion <DMIL) 

DMIL(x > y)=triie when DMV(x,y)oO AND 
Iesion(x,y)=true 



Step 4 - Measure texture of dermal melanin 
distribution (TDM) in the lesion 

TDM - SD of DMVfoy) where 



Step 5 - Calculate melanin texture ratio 
(MTR) for the lesion 

MTR = TDM/TTM 



y 

Step 6 - Threshold MTR to identify lesions 
with a high possibility of malignancy 

Malignant true if MTR > threshold 



Hgure 11 
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Figure 12 
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Step L Take spectral measurements over 
an array of points across the skin- 



. * 

Step 2. Calibrate each spectral 
measurement to eliminate effects due to 
variations in collagen concentration at 
that point, as described in GB96240O3.1 

i zzz: 

Step 3. Map each spectral BCieamjrement 
into a two dimensional colour space 
chosen so that variations to blood 
concentration in the tissue have no effect 
m tha t ai>ace. 

E 

Step 4, Using equation 1 plot a aeries of 
contours of increasing dermal melanin 
concentration st the papiltay dermis 
depth, as shown in figure 5, in that space. 



Step 5. Determine which contour is 
cfosest to the point represented by the 
calibrated spectral measurement The 
dermal roelantn concentration value of 
that contour is then used as an 
approximate value for fo(dermal melanin 
depth, dermal melanin concentration, 
epidermal melanin concentration) 



Step & Make an array of the value of this 
function corresponding each 
measurement location. 



Figure 13 
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